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We discuss the effect of the time evolution of extra dimensions on CMB anisotropies and large- 
scale structure formation. We study the impact of scalar fields in a low-energy effective description 
of a general class of brane world models on the temperature anisotropy power spectrum. We show 
that when the coupling between these scalar fields and matter evolves over cosmological timescales, 
current observations of the CMB anisotropies can constrain primordial values of the fields in a 
manner complementary to local, late-time tests of gravity. We also present the effect of these fields 
on the polarization anisotropy spectra and the growth of large-scale structure, showing that future 
CMB observations will constrain theories of the Universe involving extra dimensions even further. 
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I. INTRODUCTION 

The dominant force in the universe is gravity. According to modern cosmological theories, primordial density 
perturbations were generated in the very early universe by quantum fluctuations, whereas gravity shaped the structures 
of the universe [1]- [3]. This standard picture of structure formation by gravitational instability has gained remarkable 
support from observations of the fluctuations in the cosmic microwave background (CMB) radiation [4]. In addition, 
predictions of hierarchical clustering on scales up to 1000 Mpc in the context of cold dark matter theories (CDM) 
are in very good agreement with observations of the large scale structures in the universe. These observations may 
be used to determine the present cosmological parameters, such as the density of cold dark matter, the age of the 
universe and the equation of state of dark energy among others. 

These observations allow us not only to determine cosmological parameters, but also to test our theories of gravity 
itself. So far, no observation of a deviation from General Relativity (GR) has been reported, but, as we shall see, 
observations in the CMB and the large scale structures (LSS) provide important complementary tests to experiments 
on Earth or in the Solar System which constrain deviations from GR. The growth of perturbations predicted by a 
given theory of gravity is sensitive to the details of the theory; therefore, the study of perturbations can give important 
insights about any deviations from GR at different cosmological epochs. 

One well-known class of models where deviations from GR are predicted is that of scalar-tensor theories, in which 
the gravitational sector contains not only a tensor field but also scalar fields. The Brans-Dicke theory is an example 
of such a scalar-tensor theory (see [5] for a recent review on scalar-tensor theories). One particular feature of 
these models is that Newton's constant is no longer constant in the cosmological evolution; alternatively, in the 
so-called Einstein frame, which can be obtained from the original theory after a conformal transformation and field 
redefinitions, Newton's constant is truly constant, but the masses of particles are no longer independent of the space- 
time coordinates. Calculations of CMB anisotropies in scalar-tensor theories [6]- [8] show that the positions and 
amplitudes of the acoustic peaks usually depend on the parameters of the theory. There is no general trend on 
how the peaks are affected, because the dynamics of the fluctuations depend critically on the evolution of the scalar 
field and its coupling to matter; but the positions and amplitudes of all peaks are affected. Similarly, the slope and 
amplitude of the matter power spectrum are affected. In [6], [7] and [8] certain models of scalar-tensor theories and 
their cosmological consequences are discussed in considerable depth. In [9] and [10] models are discussed in which the 
scalar field plays the role of dark energy and couples to dark matter. In [11] it is shown that CMB anisotropies can 
probe deviations from the standard Friedmann equation. 

The studies of scalar-tensor theories and their cosmological consequences are of considerable interest, because 
theories beyond the Standard Model of particle physics imply that GR is not fundamental. Prime examples of theories 
in which deviations from GR are expected are those involving supergravity and superstrings (or their extension: M- 
theory) [12]. In particular, models which predict the existence of extra dimensions can at low energies usually be 
described as scalar-tensor theories. Typically, the low energy effective description of higher dimensional theories 
contains many light scalar fields (so-called moduli fields), which couple to matter in a manner dependent on the 
details of the higher dimensional theory. In this paper we shall discuss the scalar fields appearing in brane world 
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models and their effect on the evolution of cosmological perturbations (for recent reviews on brane worlds see, for 
instance, [13]- [15] and references therein). The brane world model we consider is of the two brane kind. The model 
is quite generic and contains the well-known Randall-Sundrum brane world model (RS I) [16] as a special case. As 
it turns out (see section II below), the model predicts the existence of two scalar fields at low energies. One of the 
scalar degrees of freedom is associated with a bulk scalar field, i.e. a scalar field which can propagate in the higher 
dimensional spacetime between the branes. The other scalar degree of freedom is related to the physical distance 
between the branes. The two scalar fields couple to matter differently and thereby affect the growth of perturbations 
in the universe each in their own way. 

In the theory to be discussed, the coupling function of one of the scalar fields is constant, its value depending on 
only one free parameter of the higher-dimensional theory. This parameter has to be chosen to be small in order for 
the theory to be consistent with local (Solar System) experiments. The other coupling function, however, depends on 
the second field itself. Barring any stabilising mechanism, this field evolves in time due to the cosmological evolution. 
Thus, even though the coupling must be small today, it could have been larger in the early universe - in particular at 
the epoch of matter-radiation decoupling. Interestingly, in our theory the matter coupling of the second field is driven 
toward small values during the matter dominated epoch, meaning that the theory can easily be made consistent with 
experiments in the solar system 1 . 

In this paper we discuss the impact of the moduli fields on the evolution of cosmological perturbations. We compare 
the predictions to the ACDM model. We furthermore assume that there is no potential energy for the scalar fields, 
so that the energy density of the scalar field is always much smaller than radiation or matter energy density. In the 
models used in this paper, the reported acceleration of the universe at low redshift is always caused by a cosmological 
constant. This ensures that the distance to the last scattering surface and the time of matter-radiation equality is 
only marginally affected by the presence of the scalar fields. Furthermore, in order to distinguish between the traces 
of the different fields we study them individually. 

The paper is organised as follows: In section II we briefly review the five dimensional setup of a brane world 
theory and discuss the effective theory in four dimensions; readers familiar with brane worlds or not interested in the 
details may skip this part. In section III we write down the action to be discussed in this paper and derive the field 
equations. We also study the cosmological background evolution of the fields. In section IV we derive the perturbation 
equations and discuss the influence of the individual fields on CMB anisotropics and the matter power spectrum, and 
we conclude in section V. This paper is addressed to readers with different backgrounds, so we have attempted to be 
comprehensive, and have included some pedagogical explanations to highlight salient points; interested readers should 
of course consult the various cited review articles for further clarification and discussion. 



II. A FIVE DIMENSIONAL THEORY AND ITS EFFECTIVE FIELD THEORY AT LOW ENERGIES 

As an introduction to the theory studied in this paper, we begin with a discussion of a higher-dimensional theory 
motivated from brane worlds. In these kind of theories, the standard model particles (and maybe also some other, ex- 
otic, form of matter) are confined on a four 2 -dimensional object (called a brane), which moves in a higher-dimensional 
space (known as the bulk). The brane itself carries some intrinsic tension as well. In our model, the bulk spacetime is 
five-dimensional. This setup is well motivated by recent developments in string theory and provides a useful starting 
point for more complicated models [18]. Also extracted from string theory is the fact that there is a second brane 
somewhere in the bulk. Confined on this second brane, which has an intrinsic tension too, can be some form of matter 
whose only direct interactions with the matter on first brane are mediated via gravity. This matter is a potential dark 
matter candidate. 

The physical spacetime stretches between the two branes, so that the branes form the boundary of the five- 
dimensional spacetime. The space between is not empty; in general, some form of matter is expected to propagate 
through the bulk: in the model under consideration, a scalar field is present. Taking this scalar field into account, the 
overall model setup is presented in figure 1. 



1 Attractor-like behaviour has been found in other scalar-tensor theories as well [17]. 
2 Three space and one time. 
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FIG. 1. The setup of the two brane model, ip denotes the bulk scalar field (with potential energy U(ip)), which can propagate 
between the branes. p\ and p2 denote the total (i.e. tension + matter) energy density confined on the individual branes. D is 
the physical distance between the branes. 

For energy densities much higher than the brane tension it has been noted [19] that the expansion of the four 
dimensional brane universe is not governed by the usual four-dimensional Einstein equations, according to which the 
expansion rate H is proportional to the matter density p; instead, it was found that the dependence of H on p is 
H oc p. This behaviour was observed in one-brane scenarios [19] but was shown to hold in two-brane scenarios as 
well (see in particular [20] and [21]). However, for densities much less than the brane tension, the usual law H oc ^fp 
can be recovered, albeit with some corrections coming from the scalar fields associated with the bulk dynamics. 

As already mentioned in the introduction, the effective field theory at low energies (i.e. at energies much less than 
the brane tension) can be written as a scalar-tensor theory with two scalar fields. The action of the theory was 
derived in [22] and reads 

12a 2 
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In this equation, 1Z is the Ricci scalar, ip and R are the two scalar fields and V is the total potential energy of the 
fields, a is a parameter of the higher-dimensional theory, which will be left free, though note that models of particle 
physics, in general, will make predictions about its value [23]. In the Einstein frame, each individual matter form on 
the branes usually couples to the fields ip and R, so that the action has the general form 

^Matter = SWtter.l (V'l > -^(V) #)flV) + ^Matter^^, B((p, R)g^v), (2) 

where ipi denotes the matter form on brane 1 and -02 the matter form on brane 2. The functions A and B in general 
depend on both fields and are not equal. 

The interpretation of (p and R is not straightforward, but the fact that there must be two scalar degrees of freedom 
at low energies can be understood as follows: first, the bulk scalar field can propagate along the brane directions, 
because it depends on all spacetime coordinates. The other degree of freedom which appears in the low energy effective 
action is related to the distance between the branes, which is a scalar function of the four spacetime coordinates. We 
note furthermore that the case where a is zero (corresponding to the second, field-dependent coupling we consider 
below) subsumes the Randall-Sundrum two brane model, where there is no bulk scalar field. 

The matter couplings are determined by the functions A and B, or, to be more precise, by their derivatives with 
respect to ip and R. For the brane world model we consider, the coupling functions are given by (see [22] for details) 

(1) = dhiA 2a 2 (2) = dlnB 25 2 

P * ~ dip l + 2a 2 ' P * ~ dip l + 2a 2 ' K) 

(i) _ 9 In A tanhi? ( 2 ) _ d\nB (tanhi?) -1 

( r = ~dR~ = 1 + 2a 2 ' 1r =~dR~ = 1 + 25 2 ' ( ' 
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The functions depend on the details of the higher dimensional theory, as they provide all information about the 
nature of the higher- dimensional space (such as its geometry or curvature) relevant at low-energy. Throughout this 
paper we shall consider coupling functions of the form above. 

There are different ways of obtaining the effective theory at low energies. In [22], the moduli space approximation 
was used, in which one integrates the heavy Kaluza-Klein like excitations out. Strictly speaking, the action above 
is only valid if the Kaluza-Klein excitations, corresponding to massive particles in four dimensions, are negligible. 
For example, massive partners of the gravitational zero mode, described by g^ u , are neglected. Another method to 
derive the low-energy effective action was used in [24]- [27], which is based on solving the five-dimensional Einstein 
equation to linear order, and then projecting onto the brane on which the standard model particles are confined. The 
theory is then written in the Jordan frame. This method is widely used in works on brane cosmology and has the 
advantage in that it provides a geometrical framework for the full five-dimensional problem 3 . However, it should be 
noted that the action with the two methods lead to the same low-energy effective theory: after a transformation to 
the Einstein frame, the action obtained by the projective approach is equivalent to the action (1). 

In practical terms, the use of the moduli space approximation limits the backwards extrapolation in time; for the 
use of the moduli space approximation to be valid, we require that the Hubble scale be much greater than the scale 
under consideration. So, H(z) w (1 + z) 2 H must be much less than fc _1 , the size of the extra dimension, so we can 
see for instance that millimetre (about 10 13 GeV _1 ) scales can be dealt with using this approximation to redshifts 
of about 10 15 (while numerical integration, of which more discussion in section IV, typically starts at redshifts of 
10 9 ). As already said, in deriving the action above, one assumes that Kaluza-Klein excitations are negligible at these 
redshifts. For example, the bulk scalar field in this approximation does not depend on the extra dimension. We 
further assume that such excitations are not produced at the epochs which we are interested in. This should be valid 
because any excitation corresponds to a heavy particle. 

Clearly this theory, described by the action in equation (1), is different from GR, and therefore gravitational 
experiments constrain parameters of the theory. In order for the theory to agree with experiments on earth and 
observations in the solar system, the parameter a has to be rather small (a < 10~ 2 ) and the field R today, should be 
small (R < 1CT 2 ), too. However, whereas the parameter a has to be chosen small from the beginning, the field R is 
dynamical and so one has to choose its initial condition such that it is small today. Given this, there is another, more 
intriguing possibility: that during the cosmological history the field is driven towards small values. The cosmological 
evolution of the system was discussed in [22], where it was shown that R is indeed driven toward small values in 
a stable fashion. Apart from local experiments, the theory is constrained by cosmological considerations. Firstly, 
nucleosynthesis bounds have to be respected, which put constraints on both a and the initial value of the field R. It is 
beyond the scope of the paper to constrain the theory by nucleosynthesis as a detailed analysis with two fields involved 
has to be performed. The second important cosmological constraint on the fields is enforced by CMB anisotropies and 
large scale structure considerations. We will in this paper develop an understanding of the how the CMB is influenced 
by the fields. 

Before we discuss perturbations in models of the form above, we comment on the frame to be used. The issue of the 
frame and its physical relevance is not new, and has been hotly discussed in the literature (see e.g. [30]). The reason 
for choosing the Einstein frame instead of the Jordan frame is that, written in the Jordan frame, one the of the scalar 
fields has a negative kinetic energy [22] (the metric in the Jordan frame is the induced metric of the first brane, where 
the standard model particles live). However, this is not a signal of an instability, because in the Einstein frame the 
kinetic term is positive. The fact that the energy density of the scalar field is not positive definite in the Jordan frame 
is not unique to field theories derived from brane worlds but appears often in scalar-tensor theories. Several authors 
have used this as an argument for using the Einstein frame as the physical frame (which is unique), and we follow 
this argument here. We would like to point out, however, that in the case of the field with non-constant coupling, 
the matter coupling of this field is very small today, so that at the present epoch the theory is indistinguishable from 
General Relativity. 

III. COSMOLOGY AND MODULI FIELDS 

As discussed in the last section, there are usually several scalar fields in the low-energy effective action, each of 
which couples to matter. The goal of this paper is to understand the effects of the individual fields on the evolution 



3 In Ref. [28] this formalism was used to calculate the CMB anisotropies in the one-brane scenario of Randall and Sundrum 
[29]. 
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of perturbations. We shall henceforth discuss them separately, dealing with only one scalar field at a time. We shall 
see in the next section that there arc differences in the predicted CMB power spectrum due to the fact that the field 
coupling function is in one case constant whereas in the other case depends on the field. We shall study two cases 
individually: the matter coupling will either be of the form (3) or (4) . 

The toy model we consider in this paper is described by the following action: 



/ dAxV ~ g & ~ 9" v (M)(9 v ct>) - V(<f>) 



+ S v {4>v,A((f>)g^) + S IV (ipiv, Bifyg^) 



(5) 



In the context of brane worlds, the index V stands here for matter on brane 1, whereas IV stands for matter on branc 
2. The functions tpv and ipiv stand for the matter fields on each individual branes. The field <j> in equation (5) takes 
the role of either if or R in the last section. In the following we shall work with reduced Planck units k = 1. 
From the action above, one can derive the following field equations: 



(6) 



n<f> = -(3 V T V - (3 IV T 



IV 



(7) 



T^ = (3 (i) T^^, with i = (V,IV). 



(8) 



In these equations, T stands for the trace of the energy-momentum tensor for each individual matter. The functions 
(3 are given by 



d\nA 



02 = 



d\nB 



(9) 



In the first case under consideration, the functions (3 do not depend on the scalar field (as it is the case with one 
of the scalar field in section II) and choose (3 to be the form 4 



Case I : A = B cx cxp(-2a 2 0/(l + 2a 2 )) 



(10) 



We have retained the main features of the brane model, i.e. [3\^ is negative and bounded. In the second case 
(corresponding to the field R in the last section), the functions (3 depend on the field and we have 



Case II : A cx cosh(</>), B cx sinh(</>) 



/?i = tanh(^) = 1/02. 



(11) 



As before we have kept the same functional dependence as in the brane models. 

For a homogeneous and isotropic universe with flat spatial sections, the line element reads 



The field equations read 



ds 2 = a 2 (r) (-dr 2 + 5 i jdx i dx j ) 



(12) 



H2 = l a2 (pv+Piv + ^ 2 + V(cf>) 



(13) 



• + 2H(j> + a : 



dV 



-f3v(pv ~ 3p v )a 2 - f3iv(piv - 3piv)a 2 



(14) 



P(i) + 3i?(p(j) +P(i)) = (3(i) (p(i) - 3f> (i) ) (j). 



(15) 



'Note that a in this expression is different from a in Section 2, because of the field redefinition. 
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In these equations, the dot represents a derivative with respect to the conformal time t and H = a/ a. Note that in 
equations (7) and (14) we must sum over all matter forms. In this paper we will assume that = 0. This would 
correspond to no matter on the second brane in the theory of the last section. 

The theory above differs from General Relativity and is therefore constrained by observations. Firstly, local ex- 
periments (i.e. on earth and in the Solar System) constrain the post-Newtonian parameters 7pn and /3pn given 

by 



B 2 



A/B 
A/B 



7pn-1 = -2— (16) 



_ 1 Pa/b d(3 A/B 
/?PN - 1 -2(l + /3i /B ) 2 # ' (17) 

where /?a/b stand for the functions (B\ for the cases I or II, respectively. Current constraints give [31] 

|7PN-1| < 2 • 10- 3 , |/3pn-1| <6-HT 4 . (18) 

Nucleosynthesis constraints limit the effective number of degrees of freedom for relativistic particles at this epoch. In 
general the constraints have to be worked out in detail, but there is a simple way to get a rough limit on the theory 
from nucleosynthesis [32]. Because the energy conservation equation (15) implies that Pmo 3 constant in general, 
the expansion rate at nucleosynthesis is different from its value in General Relativity. This leads to [32] 



< 1.2. (19) 



Although we have derived the general equations, from now on we shall set V{4>) = A, because we wish to compare 
the results for the perturbations with the ACDM model. Before we turn our attention to cosmological perturbations 
and their evolution, however, we briefly study the evolution of the background in the matter dominated era for the 
case A and B separately. 

Case I: The case for constant coupling parameter (3 

The model in this case has similarities to the model discussed in [33], although we do not have an exponential 
potential for the field, and so the critical points found in [33] do not apply here. In the radiation dominated epoch, 
the field is almost constant, because the terms on the left-hand side on the Klein-Gordon equation (14) dominate. 
The constraints (18) on the post-Newtonian parameter lead to a < 0.1. We can easily find the solutions to the 
background equations in the matter dominated epoch: 

o(r) = ^— j , with x = 2 - 8a 4 

= 0o + 4a 2 In (^j . (20) 

Therefore, the universe expands a little slower than expected, although the correction is very small. An important 
point has to be made here, namely that matter does not scale like p ma ttcr oc a~ 3 , but rather like 



Pmattor OC a 3 exp 



J • (21) 



For the theory under consideration here, this implies that, for a given matter density today, the matter density in 
the past was greater than it would have been under standard cosmology. 

Case II: The case for non-constant coupling parameter (3 

The low-energy effective theory now corresponds to the first model by Randall and Sundrum [16]. As we will see, 
there is an attractor mechanism at work, which drives the field towards small values and thereby generating small 
couplings between <f> and matter, as it is dictated by observations on Earth and in the Solar System. 

The equations are difficult to handle analytically. Assuming that <j) is small and slowly varying at the onset of 
matter domination, one can find that 
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(22) 



and that p ma ttor oc a 3 . The field <j)(r) now goes like 



cos(Vf ln(r/r )) sin (Jf In (r/r )) 



(r/ro) 3 / 2 



(r/ro) 3 / 2 



(23) 



Thus, apart from oscillating behaviour, the field tfi in this case decays during matter domination. 

Although the solutions above are based on a crude approximation, this qualitative behaviour is found even if the 
field value was initially larger [22]. For large values of <f>, however, equation (21) is a better approximation. 

Thus, although it is easy to fulfill the constraints (18) due to the cosmological attractor, the nucleosynthesis 
constraints give a limit on <j>(r = r nuc i) < 0.4. 

IV. COSMOLOGICAL PERTURBATIONS AND MODULI FIELDS 

In the following we discuss the evolution of perturbations in the presence of moduli fields. This section is organized 
as follows: first we will write down the perturbation equations. Then we will solve them (approximately) in the matter 
dominated era, in order to gain some understanding about the effects of the coupling between the field and matter. 
Finally we discuss the solutions of the numerical computation of the spectrum of anisotropics in the CMB and the 
matter power spectrum. 



A. The Perturbation Equations in the Synchronous Gauge 

We work in the synchronous gauge (using the notations as in [34]); i.e. the perturbed line element has the form 

ds 2 = a 2 (r) (-dr 2 + {5ij + h^dx'dx 1 ) . (24) 
After Fourier transformation, the perturbed Einstein equations read: 

(25) 
(26) 

where h and r\ are defined by (k is the unit vector in direction of k) 



2k 2 r]- -h = a 2 ST° Q 
a 

h + 2-h - 2k 2 r] = -a 2 5T\ 



hij 



I 



d 3 ke ikx 



kj • kj/i(k, t) + ( kj • kj - -Sij ) 6»?(k, r) 



(27) 



^From the energy-momentum conservation equation one obtains a set of two equations for the evolution of the 
density contrast and the divergence of the velocity field. Defining 



(p + p)0 = ik 3 ST° 3 

(p + p)a = -(k- k ^) (T t] ^n/3) 



(28) 
(29) 



and writing pi = c 2 ^ = dp^/dp^ and 9 = ik 2 Vi we find from the evolution equation for the density contrast 

S(i) = 5p(i)/p(i) 



<5 W =-(l + w {i) ) 



§| (1 " Mo) 



Pi (l - 3w w ) (8<j>) - 30/3, (c 2 s{t) - 8 {{) . 



(30) 
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The divergence of the velocity field 9^ = ■ of each individual fluid evolves according to 



%) = -- (1 - 3w (t) ) 6 {i) - j^-O^ + —^k 2 S - k 2 a (€) 



s(i) 



1 - 3w 



(0 1.2. 



w 



« 



+ w {i) 1 + w (l) 

A(l - 3tu (i ))^ (i ). 



(31) 



It is important to note that these equations reduce to those of ordinary General Relativity, if the fluid under 
consideration has an radiation-like equation of state, i.e. if wu\ — 1/3. This means that the equations governing the 
dynamics of photons and relativistic neutrinos are those found in General Relativity. Important changes occur in the 
equations for CDM and baryons: 

Cold Dark Matter: 
In this case, we have w = = c s and vanishing anisotropic stress. Thus 



(32) 



C = —6 c + p c k 2 54>~p c <t>0 c 
a 



(33) 



Baryons: 

In this case, we have w = = c s and vanishing anisotropic stress (a very good approximation after neutrino decoupling 
[34]). As in [34], we do not neglect the quantity c 2 k 2 5, which becomes important on small length scales and thus large 
k. We also have to add the terms describing the momentum transfer between photons and baryons due to Thomson 
scattering. Therefore, wc have 



(34) 



ti,, - ~9 b + c 2 s k 2 S b + ^an e o T (0 1 



3^6 



9 b ) + Pbk 2 



(35) 



where n e is the number density of free electrons and &t the Thomson cross section. 

The only new terms compared to General Relativity are the ones proportional to (3 and its derivative. Note that 
we have to add an evolution equation for 9 C : while we can still set ti c = initially (in order to fix the synchronous 
gauge completely), in General Relativity, where there is no extra field, it stays zero; however, fluctuations in <f> are a 
source and therefore, in the theory we consider, it evolves according to the above equation. 

Finally, the equation for the fluctuations in <j> evolve according to 



( M+2^)+(fc 2 + a 2 ^ 



5<p + -h<j> 



These equations are quite general. Again, we are assuming V = A = const. 



B. Solutions to the Perturbation Equations in the Matter Dominated Era 



(36) 



Before we give the results for the CMB anisotropies and matter power spectrum, we try to gain some analytical 
insight from the perturbation equations in the matter dominated epoch on sub-horizon scales. We neglect radiation 
and the baryons, taking only into account CDM and the scalar field in this regime. Furthermore, we consider only 
the case for the field with constant coupling, which is tractable analytically. 

Taking the derivative of equation (32) and making use of equations (25), (26), (33) and the Friedmann equation, 
one gets 
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S c + H5 C - ^H 2 5 c = 2a 2 k 2 S<j) + 20(50) - 2a 2 (H(8<t>) + (<ty)) - 2a 2 ^6 c . (37) 
Using the perturbed Klein-Gordon equation, ignoring oscillations in 5(j) and inserting (32) one obtains 

S c + (H - 2a 2 <j>)6 c - ^H 2 (l + 8a 4 ) 6 C = 0. (38) 
Assuming a cx t x , using the background solution equation (20) and making the Ansatz S c oc r m , one gets up to 0(a 4 ) 

m = 2 + %-a 4 (39) 

for the growing mode. Therefore, compared to the ACDM model, due to the non-vanishing coupling perturbations 
in CDM grow faster than the scale factor and the transfer function is on small scales bent towards larger values. We 
have confirmed this numerically. 



C. CMB Anisotropies and the Matter Power Spectrum 

In the following we will discuss the spectrum of anisotropies in the CMB and the matter power spectrum. It is 
customary to expand the temperature anisotropy measured in a given direction n in spherical harmonics according to 

AT 

— (n) = ]Ta ;m YUn). (40) 

Ira 

We assume that the density fluctuations are Gaussian random variables, as predicted by the inflationary scenario. 
The statistics of the fluctuations are completely specified by the power spectrum Cj, defined as 

(ai rn 

a l'm') = Cl (41) 

Let us briefly summarize how the C/s depend on the perturbation variables and how they are calculated (see [35], 
[36] and [3] for excellent detailed reviews on the physics of CMB anisotropies, [37]- [40] are shorter reviews). The 
reader familiar with CMB anisotropies can skip this part. To illustrate the essential steps involved in calculating the 
anisotropy power spectrum, we neglect the polarization states of the photons. However, when we present the results 
of the calculations, the polarization states have been taken into account. 

At very early times the universe was radiation dominated and the temperature of the universe was so high, that 
photons and baryons formed a plasma. At these high temperatures, photons and baryons are tightly coupled. On 
scales less than the horizon the fluid undergoes acoustic oscillations due to the interplay of gravity, which tries to 
compress overdense regions, and the photon pressure, which, when high enough, acts against the force of gravity. 
When the temperature dropped below 4000 K, electrons and nuclei were able to combine to atoms (mainly hydrogen 
and helium) and the mean free path of the photons became larger than the Hubble horizon. Around that time, the 
photons completely decoupled from matter. On their path, photons may experience energy loss due to time-changing 
gravitational fields (the so-called integrated Sachs- Wolfe effect, or ISW). Additionally, photons may interact with 
reionized gas (reionization might happen due to early star formation). 

In order to follow the perturbations in the photons correctly, one has to go beyond the fluid description and study 
the evolution of the photon distribution function / 7 in phase space, i.e. one has to solve the Boltzmann equation. In 
general, / 7 is a function of the comoving position x, the photon momentum p and conformal time r. To first order we 
write / 7 = /o + Sf 7 , where /o is the unperturbed photon distribution. The brightness function Aj(x, n, r) is defined 
as 

<y/ 7 (x,p,r)=(^)A T (x,n,r), (42) 

where T is the unperturbed temperature of the photon gas (i.e. of the background) and n is a unit vector in the 
direction of the photon momentum p. One can now formulate a perturbed version of the Boltzmann equation, 
describing the evolution of At- However, it is more useful to perform a Fourier transform of At(x, n, r) and expand 
the Fourier coefficients A^(k, n, r) in Legendre polynomials (k is a unit vector in direction of k): 

oo 

A T (k, n, r) = ^(-i)'(2/ + l)A T/ (fc, r)P,(k ■ n). (43) 
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It can be shown that the density contrast in the photon energy density, 8 1 , the divergence of the photon velocity field, 
8 7 and the photon anisotropic stress, tr 7 , are related to At; by 

3 1 

<5 7 = A T0 , 6 1 = -kA TU cr 7 = -A T2 . (44) 

In terms of the moments At;, the Boltzmann equation becomes a hierarchy of equations. Using equation (44), the 
evolution equations are given by 



8 7 = k 2 — °7^ an e cfT (Ob — # 7 ) , (46) 

8 3 4-8 9 

2a 7 = — 9 1 - -kA T3 + —h + -f) - -em e <7 T cr 7 (47) 

k 

A « = ^7 [/A T( ,_i) - (Z + l)A T( ,+i)] -an e a T A T i, (I > 3) (48) 

Furthermore, it can be shown that the quantity Ci, defined in equation (41), is related to the brightness function At; 
by 

G=4^|*|A T ;(fc,r )| 2 . (49) 

It is straightforward to include polarization. Due to Thomson scattering, photons are linearly polarized in the plane 
perpendicular to n. The collision terms in the Boltzmann equation depend on the polarization state. In order to 
follow the evolution in phase space, one defines a total distribution function / 7 , which includes all polarization states 
and another distribution function, g 7 , which is defined as the difference of the distribution functions for the individual 
polarized photons. Defining the temperature brightness function At and the polarization brightness function Ap 
similarly to the above procedure, one can show that At obeys the equation (45)-(48), but equation (47) has to be 
changed to 

8.3,, 4-8.9 1 , A AX 

2cr 7 = — 7 - -fcA T 3 + j^h + -7] - -an e a T o-j + —an e a T (A P0 - A P2 ) , (50) 

because the collision term in the Boltzmann equation depends on the polarization state. Additionally, one obtains a 
hierarchy for the functions Ap;: 



Ap; - ^j-j- [/A P(i _ 1) -(1 + l)A P(i+1) ] - an e a T A Pl + ^ (A T2 + A P0 + A P2 ) (s M + ^ 



(51) 



where Su in the last term stands for the Kronecker delta. 

The CMB anisotropy power spectrum C; is still given by equation (49). Thus, in order to calculate the CMB 
power spectrum C; (as well as the polarization power spectrum and the temperature/polarization cross correlation 
spectrum), one has to solve a large set of differential equations (the equations for baryons, CDM and the scalar field 
presented in section IV A and the hierarchy above). In [41] solutions to the Boltzmann equation are found. Here we 
use an existing code [42] and make the appropriate changes to the evolution of background and perturbation variables 
to incorporate the above theory 



1. Constant Coupling 

In order to calculate the CMB power spectrum, we make changes to an existing code [42], whose method is based 
on [43]. As described in section IV A we have to make changes in the equations for CDM and baryons and have to 
include the modified Klein-Gordon equation for the scalar field. We assume a scale-invariant initial power spectrum 
and adiabatic initial conditions. 

What are the expected effects of the non-vanishing coupling? There are two effects, which come in when a differs 
from zero. Firstly, as mentioned in section III, the dependence of the matter density (both for baryons and CDM) on 
the scale factor is no longer given by p ma ttcr oc a~ 3 , but rather by equation (21). Therefore, relative to the ACDM 



10 



model the density of matter (and in particular of baryons) is larger at decoupling (we remind the reader that the 
models are normalized such that their matter density is equal at the present epoch). Hence, the amplitudes of the 
peaks are larger, since a larger baryon content implies larger amplitudes of the peaks. 

Additionally, there are changes to the sound speed c 2 = Ay/(Ab + Ay)- The sound speed decreases with increasing 
a. This effects the peak separation, which is proportional to [35] 



61 = 



X(T0 ~ Tdec) 

D x 



(52) 



where x( T o — <dcc) is the angular diameter distance to the surface of last scattering and D s = f Tdoc Cs dr is the sound 
horizon. Aside from having an effect on c s , the presence of the scalar field does also affect the distance to the last 
scattering surface. Let us estimate the position of the first peak, given by Z = 7r x( T o — T& ec )/D s . It is simple to 
calculate x(t~o - Tdcc): 



/"TO 

x(t - Tdcc) = / At ■. 



1 



H l + C 

where C = 2a 4 /(l — 8a 2 )(l + 2a 2 ). Putting c 2 w 1/3 and neglecting the radiation dominated epoch, one finds 



1 1 



Therefore, the position of the peak is given by (note that we have a(ro) = 1) 

7 _ 1 - C / 2 

to - toa d cc > 



(53) 



(54) 



(55) 



where Zo ~ 200 is the position of the peak for the flat model with no coupling. Therefore, the peak moves towards 
larger values (i.e. to the right in the spectrum). 

To reach this conclusion we have assumed that c s w 1/3, which is a crude assumption. Note, however, that if one 
takes the baryons into account the c s is smaller and hence, the sound horizon D s is smaller. Thus, Zo is if anything 
somewhat larger than the value given in equation (55). 




FIG. 2. The temperature anisotropy power spectrum, 1(1 + 1)Ci/2tt, for the constant coupling case: the values in the legend 
are the values of a. 



The second important change when the coupling is non-vanishing is that there is a larger contribution from the 
integrated Sachs-Wolfe effect. Even in the matter dominated era the gravitational potential changes with time, 
because the matter density contrast no longer grows like S(t) oc o(t), as discussed in the last section. As the time 
evolution of the scale factor is altered in the presence of the scalar field, the duration of the recombination era is 
affected, too. 
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In figure 2 we plot the results for the temperature anisotropy power spectra. We assume an initial power spectrum 
with equal power on all scales. It can clearly be seen that the both the amplitude and the separation of the peaks is 
affected with increasing coupling. Also, on large scales (low I) one can see the additional power due to the integrated 
Sachs- Wolfe effect. 

In order to compare the theory to observations, we normalize the curves appropriately. Usually this done at large 
angular scales, using the COBE data as reference. The curves normalized to COBE are plotted in figure 3. Clearly, 
the normalized curves show less power on small scales because of the enhanced integrated Sachs- Wolfe effect in 
models with non-vanishing coupling. On larger scales (small multipole number I) one clearly sees a boost of power. 
Additionally, the normalized matter power spectra for models with non-vanishing coupling are below the ACDM 
curve. 

For completeness we plot the E-type polarization spectra and the cross correlation between temperature and E-type 
polarization in figure 4. One can clearly see that in both spectra on small scales the power is suppressed when a is 
increased. 
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FIG. 3. COBE-normalized temperature anisotropy 1(1 + 1)C\/2-k (left panel) and matter power spectrum (right panel) for 
the case of constant coupling. On small scales the COBE-normalized spectra are below the predictions for vanishing coupling 
due to the enhanced ISW. 




10000 



FIG. 4. COBE-normalized E-mode polarization anisotropy (left panel) and TE cross-correlation (right panel) for the case of 
constant coupling. 
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FIG. 5. CMB anisotropy spectra (left panel) for the constant-coupling case, with a taking the extreme value permitted by 
nucleosynthesis, with a reference curve from the ACDM model; the right panel shows the relative discrepancy. Although the 
differences are around 3 percent at high multipoles, a different choice of normalization - normalizing to COBE at I = 800, say 
- would move all the observed discrepancies into the low-multipole region, where cosmic variance dominates. 

In figure 5 we plot the anisotropy spectra for the ACDM model and the model with the maximal value of a allowed 
by nucleosynthesis (a = 0.1). One can see that the COBE normalized curves differ from each other by a few percent 
on small angular scales. 

2. Non-constant Coupling 

In this case, where the field 4> corresponds to the R field in equation 1, the coupling function now decays. But also 
here, the baryon density is larger at decoupling compared to the ACDM, which leaves a similar effect on the peak 
position and amplitudes as in the case of constant coupling. 




FIG. 6. The temperature anisotropy power spectrum, 1(1 + l)C'i/2ir, for the field-dependent coupling case: the values in the 
legend are for the initial values of the scalar field (j>. 



However, well inside the matter dominated epoch, the coupling becomes very small, so that the density contrast 
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grows essentially like a(r) from that time on. Therefore, the ISW effect is not as pronounced as in the case of constant 
coupling. Only if (j> was initially very large will the field have not evolved to small values by today; consequently, for 
the field-dependent coupling case, only for large initial field values will the effect on the ISW become more apparent 
(see figure 6). 
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FIG. 7. COBE-normalized temperature anisotropy 1(1 + l)G'i/2ir (left panel) and matter power spectrum (right panel) for 
the case of field-dependent coupling. Similar to the case in figure 3, on small scales the COBE-normalized spectra are below 
the predictions for vanishing coupling. 
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FIG. 8. COBE-normalized E-mode polarization anisotropy (left panel) and TE cross-correlation (right panel) for the case of 
field-dependent coupling. 



The normalized anisotropy and matter power spectra are shown in figure 7, whereas the normalized E-mode 
polarization anisotropy and the TE-cross-correlation are shown in figure 8. Increasing the initial value of the field 
implies that the power in both the E-polarization spectrum and the cross correlation spectrum between temperature 
and E-polarization is suppressed on small scales. 

The initial field values in figures (6) and (7) have been used without taking into account constraints from nucleosyn- 
thesis (19), which implies that nuc < 0.4. In figure 9 we plot the cases for 0; n it = 0.4 and the ACDM curve. One can 
see that there are differences in the predictions for the temperature anisotropy as large as 6 percent at small angular 
scales and 3 percent on degree scale. This is potentially distinguishable with WMAP and the Planck Surveyor. 
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FIG. 9. CMB anisotropy spectra (left panel) for the case of an initial value (0i n i = 0.4) allowed by nucleosynthesis. The 
right panel shows the relative error between the ACDM and the model with 0i n i = 0.4. For / ~ 100 the difference between the 
model is as large as 3 percent and gets larger at higher I (smaller angular scales). 

We would like to point out that constraints on <p become stronger, once the CMB is included. In all cases plotted 
in the figures, the field fulfills the constraints today. However, a large field value in the early universe modifies the 
predictions for the temperature and polarization power spectra (as well as the matter power spectra). Hence, the 
CMB provides vital complementary information to the local experiments. 

V. CONCLUSION 

In this paper we have investigated the implications of time-varying extra dimensions for CMB anisotropies and 
large scale structure formation. In theories with extra dimensions, certain scalar fields, so called 'moduli fields', appear 
in the low energy effective theory coupling directly to matter fields. We have investigated in particular the couplings 
appearing in a five-dimensional model based on the brane world idea. The model contains two branes which are the 
boundaries of the higher-dimensional spacetime, with, in general, four-dimensional matter fields confined on these 
branes. Additionally, we allowed for a scalar field propagating in the extra-dimensional spacetime. At low energies 
the dynamics of this system is described by a bi-scalar-tensor theory. One of the scalar fields has a constant (i.e. 
field- and coordinate-independent) coupling to the matter fields, whereas in the case of the second field the coupling 
function is field dependent. 

We have investigated the effects of the individual fields on the CMB anisotropies, the cross-correlation between 
temperature and polarization of the CMB and the matter power spectrum. Due to the coupling, both the amplitude 
and the positions of the peaks are affected. The spectra normalized to COBE lie below the spectra of the reference 
ACDM model on small scales. 

The results presented in section IV imply that, for the brane world model presented in section II, potentially there 
is a measurable effect of both fields on CMB. If it turns out that the second field needs to be already quite small 
during the radiation dominated epoch, then some mechanism has to drive it to these small values in the very early 
universe. This can happen, for example, if there is some matter on the second brane. In this case, the field R is 
driven even faster to zero, so that it might be small at decoupling. Note, however, that in this case the evolution of 
perturbations are affected in the radiation dominated epoch. 

Clearly, the CMB anisotropies as well as the polarization provide complementary information about time-evolution 
of the extra dimensions. Together with other tests, such as gravity experiments in the solar system and/or on earth 
and tests of time-varying gauge couplings and primordial nucleosynthesis, these cosmological considerations strongly 
constrain any model involving extra dimensions. For example in the case of the model considered here, the field with 
constant coupling is already constrained by local experiments such that its effects on the CMB are negligible. For 
the second field, however, we have seen that although its effect is constrained to be minimal today, it can have a big 
impact on the CMB anisotropies. 
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There is another important point not addressed in this paper: we have not discussed the time evolution of Kaluza- 
Klcin excitations and their influence on the CMB. They might still be present in the radiation dominated epoch, 
thereby affecting the evolution of perturbations at least initially. Their presence could imply, for example, that the 
initial perturbations arc not purely adiabatic. The details would depend on the mechanism which produces the 
primordial fluctuations, such as inflation or the cyclic model, see e.g. [44] and [45]. To address this important issue 
one has to go beyond the moduli space approximation (and in fact beyond the lowest order in the method used in 
[24], [25] and [27]). 

In this paper we have investigated the effects of moduli fields on the CMB motivated from a certain class of brane 
world theories. There are other brane world models, in which the theory at low energies cannot be described by an 
action of the form (1). For example, in the models of [46]- [48], gravity becomes five-dimensional at large distances. 
We expect that some of the conclusions drawn in this paper do not hold for such models. Nevertheless, the CMB 
gives useful constraints on these models, too [49]. 
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